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Abstract 

A finite one-dimensional microscopic model of a superfulid is presented. The model consists of 
interacting Bose particles with an additional impurity particle confined to a ring. Both semiclassical 
and exact quantum calculations reveal dissipationless motion of impurity with increased effective 
mass due to its interaction with the excitations of Bose fluid. It is shown that both the excitation 
spectrum of Bose fluid and the excitation spectrum of impurity can be analyzed using the structure 
of the ground state of the system. 
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I. INTRODUCTION 



Despite its maturity, the theory of Bose-Einstein condensates (BEC) and superfluidity 
remains an active field of research. New theoretical developments are motivated by a variety 
of modern experiments. Bose Einstein Condensates (BEC) of ultra-cold atoms in magnetic 
traps [H E] is currently an active experimental field. Another set of recent experiments stud- 
ies microscopic superfluidity of liquid helium by spectroscopic measurements on molecules 
imbedded in superfluid Helium droplets [31 HJ E] . It is generally observed that a microscopic 
impurity interacting with Bose liquid/gas behaves as a free particle with an effective mass 
that is greater than its original mass. Such dissipationless quantum motion is observed both 
for translational motion of particles in a superfluid environment as well as for rotations of 
molecules in superfluid helium droplets. The development of quantum microscopic theory of 
Bose fluid (BF) / impurity system has been the theoretical challenge. What is the collective 
wavefunction of this system and what is the nature of the effective mass? What happens with 
finite number of Bose particles and how is the macroscopic limit of superfluidity achieved? 
A number of theoretical works have been devoted to answering those questions. 

The system that has been accessible to analytical theory is a dilute, weakly interacting 
BEC (see, for example, [6] and references therein). In the seminal work of Bogoliubov [7j 
the excitation spectrum and the nature of excitations of dilute BEC was uncovered. This 
work set the stage for most of the further theoretical studies of this system. A number of 
authors used perturbation theory to consider the interaction of an impurity particle with 
BEC 0, El Uni EH [12] ■ It was concluded that the Landau criterion holds for the quantum 
motion of a particle in the BEC, i. e. the motion of a particle is free up to a critical momentum. 
At low momenta the particle has the spectrum that resembles that of a free particle, E(k) ~ 
k 2 , thus leading to an effective mass approximation. The effective mass was calculated by 
several authors within the Golden Rule limit of the particle/BEC interaction. In our recent 
paper we continued this work by developing a formal perturbation expansion of this system 
based on Coherent State Path Integral formulation of the particle/BEC dynamics [12]. We 
presented the diagrammatic representation of perturbation expansion that allows to evaluate 
particle properties in dilute BEC up to an arbitrary order of perturbation theory. Besides 
the perturbation theory another theoretical approach that gained popularity is the solution 
of the Gross-Pitaevskii equations for the BEC [T3] . Gross-Pitaevskii can be viewed as the 
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mean field equations of motion for the Bose field. It has been shown that many of the results 
obtained from the Gross-Pitaevskii equations, including dissipation and the effective mass 
of impurity particles, are the same as the ones obtained from perturbation theory [9]. 

Despite its important physical insight, the perturbation theory often cannot be applied 
to calculate observables for realistic systems such as liquid helium. An approach taken by 
several theoretical groups has been the computational imaginary time path integral tech- 
niques |14j . The most recent work has been aimed at the studies of molecular rotations in 
superfluid helium [Eft M> EH EE1 ESI 1201 121] • Whaley and coworkers PU EH HH1 ESI [20] 
has calculated the properties of helium droplets with a variety of molecules. Roy et al. [21] 
has considered the limit of very small helium clusters. In a number of cases the effective 
moments of inertia of molecules were computed and were shown to be in agreement with 
experimental values. However, being a powerful computation tool, this work is necessar- 
ily limited to the calculation of statistical rather than dynamical properties of the system 
and does not give a direct answer to the intriguing question of how does the motion of the 
molecule and motion of a superfluid uncouple. 

A microscopic theory that has been successful in predicting the excitation spectrum of the 
bulk superfluid helium originated from the work of Bijl and Feynman [221 123]. In this work 
the trial wavefunction of elementary excitation is built using the unknown wavefunction of 
the ground state by the use of the Fourier component of the density operator, fa = p^o- 
It is shown that the variational energy of such state is expressed through the structural 
properties of the ground state and is given by E(k) = k 2 /[2mS(k)}, where S(k) is the 
superfluid structure factor. The Bijl-Feynman spectrum has qualitatively correct shape; 
however, quantitatively it overestimates the excitation energy in the most relevant roton 
minimum spectral region. Subsequently the theory has been extended to include the multiple 
Feynman excitations, an approach generally known as a correlated basis functions (CBF) 
theory [22]- Several realizations of this theory lead to the excitation spectrum that agrees 
with the experimental one. Two type of extensions of the CBF have been applied to calculate 
the excitation spectrum of the 3 He impurity in the superfluid 4 He. One approach is based 
on the variational method [24] , while the other is based on the perturbation series [23 EE] . 
Both calculations lead to results that are in quantitative agreement with experiment. 

In this work we present a computational study of quantum dynamics of a very simple 
model system. The developed model retains many important features of the Bose superfluid, 
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it is not limited to the low density and weak interactions, yet it is small enough that its 
quantum dynamics can be computed and analyzed exactly. The calculations presented below 
clearly show the effects of microscopic superfluidity in a finite, one-dimensional system. We 
analyze the results using the analytical theory developed for BEC and superfluid helium. The 
work provides an extensive illustration of these theories and allows to clearly understand the 
limits of their applicability. In the next Section the Hamiltonian of the model is presented. 
In Section III the semi-classical equations of motion for the model are solved. It is shown 
that the model exhibits microscopic superfluidity and the nature of effective mass is revealed. 
In Section IV the exact numerical solution for a pure BF is presented. We find the excitation 
spectrum of this system and compare it to the Bogoliubov's spectrum. In Section V the 
free-particle like excitation spectrum of impurity is obtained. Finally, in Section VI we show 
that the approach based on the Bijl-Feynman theory captures the essential physics of the 
impurity excitation spectrum and provides an accurate method for calculating the effective 
mass of impurity based on the structural properties of the ground state. 

II. MODEL 

The main object in our model is one-dimensional BF. In order to make it accessible for 
numerical calculations the system is made finite by confining the BF to a ring of radius R 
(enforcing periodic boundary conditions). The system schematically represented in Fig. 1. 
The Bose particles represented by empty circles interact by pair potential U — Xj), where 
Xj and X j £1X6 the coordinates of the particles on a ring, the difference X <i X j IS the shortest 
distance between particles. An impurity which interacts with BF via the pair potential 
V(xi — x), where x is the coordinate of an impurity and Xi is the coordinate of ith Bose par- 
ticle, can be added to this system. For the sake of computational simplicity the interaction 
potentials are taken to be Gaussians of unit widths, U(x) = ae~ x2 , and V(x) = (3e~ x2 . 

Given such toy model, its exact Hamiltonian can be written in the secondary quantization 
form. The BF Hamiltonian adapted to one dimension is 

k 2 1 

**U = Yl Q^fi^k + 2L ^ Ul "t+A-i^n , (1) 

k k\,k2,l 

where U\ = f U(x)e llx dx is the Fourier transform of the pairwise interaction potential, m is 
the mass of Bose particles, and L = 2nR is the length of the ring. The summations over 
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momenta k are carried out over all possible k = 0, ±1/R, ±2/R, ... . Throughout this work 
we use a simple unit system in which ft = 1, all other variables are dimensionless and are 
of the order of unity. The total number of Bose particles, N is an input in our calculations. 
The majority of the calculations are carried out for a liquid like system, where the average 
distance between particles is similar to the range of the pair potential, N/L ~ 1. Since 
the interaction potential is smooth we do not expect large momenta to play an important 
role. It is verified in all our calculations that results are well converged when the momentum 
space is truncated by —N/R < k < N/R. 

We will consider this model both with and without impurity. An impurity Hamiltonian 
is that of a free particle interacting with the Bose particles and can be written in coordinate 
representation for the impurity, 

k,l 

For the quantum calculations it is more convenient to express this Hamiltonian in momentum 
representation, 

Himp = ^2^\p){p\ + j^2v l \p + l)(p\ at_ t a k , (3) 

V p,k,l 

were \k) = [1/ L) l / 2 e lkx is the momentum wavefunction of impurity, M is the mass of the 
impurity and Vi is the Fourier transform of the impurity/Bose particle pair potential. 

Another form of the system Hamiltonian that proves to be useful is its expression through 
the Fourier components of the density operator, 

3 I 

The BF Hamiltonian Q can then be written as 

k 2 1 

Hbf = Yl ?^fi h k + ^lJ2 U k Pk P-k , (5) 

k k 



and the impurity Hamiltonian (pj becomes 



1 d 2 1 

H imv = T^r + - VT4 e ikx h k . (6) 

imp „ , , p. 2 j £^ H k \ J 



2Mdx 2 L 

k 
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III. CLASSICAL DYNAMIC, DEMONSTRATION OF SUPERFLUIDITY 



Even though full quantum calculations are obtained in the following sections, it is in- 
structional to examine the classical dynamics of this system first. We notice that the Bose 
creation/annihilation operators of the system are equivalent to a set of oscillator degrees of 
freedom, i.e. H = |^Sfc«fc is the hamiltonian of a harmonic oscillator with the frequency 
k 2 /2m. The dynamics of such degrees of freedom can be considered classically. The in- 
teraction term presents a complex unharmonic coupling between the oscillator degrees of 
freedom. The formal use of Hamiltonian dynamics for H^f + H imp results in the following 
equations of motion for the Bose liquid: 

* = &• (7) 

dH({a* k ,a k },x,p) 

P dx ' { } 

= dH({at,a k },x,p) 

_ . dH({a* k ,a k },x,p) 

where the function H({a* kl a k }, x,p) can be obtained by replacing operators a k and a k by 
complex variables a* k and a k respectively, i.e. 

k 2 1 p 2 

H({a* k ,a k },x,p) = ^ a *k a k +2l^2 Uia* ki+l a* k2 _ iakl a k2 + — + ^2v l e dx a* k _ l a k .(11) 

k ki,k2,l k,l 

The variables a and a* can be related to the classical momenta and coordinates via a k = 
~^(x k + ip k ) and a* k = ^(x k — ip k ). The use of complex variables is convenient for such 
a Hamiltonian. The equations of motion, Eqs. rt7 10), are equivalent to the Hamilton's 



equations given that Eq. (10) is the complex conjugate of Eq. (|9j). Notice that in spite 
of the word "classical" this dynamics is not written in coordinates of particles. With the 
exception of the impurity particle the dynamics is carried out in the space of occupation 
numbers. Thus, the total number of particles is a variable. However, total momentum 
and total number of particles can be shown to be rigorously conserved. These equations of 
motion are in fact equivalent to the Gross-Pitaevskii equations [13J. 

An important drawback of this theory is that the energy minimum does not depend on 
the interaction strength and is achieved by placing all of the Bose particles in the k = state. 



An analysis of the Hamiltonian Eq. (11) shows that Hf,f(a k , a* k ) has the global minimum at 



cio = a>Q = vN, and a k = a* k = for k ^ 0. The latter was verified by performing an 
extensive numerical search. As a result such theory deals with rather unphysical state of 
the BF and the classical results may deviate significantly from quantum calculations. 

Fig. 2 shows the motion of the particle for the different parameters of the Hamiltonian 



Eq. (11) with R = 3.0. The initial conditions for these trajectories are taken to be x 



0, p = Po, do = a-o = VN, N = 19 ~ L. Of course, without any interactions the solution 
is a straight line x(t) = pot/M, shown by a dashed line in both plots as a reference. If the 
interaction between impurity and Bose gas is turned "on" but the interaction part of the 
bose Hamiltonian is zero (the case of ideal Bose gas with pair potential parameters a = 0, 
/3 — 0.5) the scattering of impurity is clearly observed in the solution as shown in the top 
plot. The scattering is rather complicated and despite the ordered initial conditions results 
in a diffusion like motion after t <~ 5. The most important result is shown in the bottom plot. 
The Hamiltonian now includes all of the interaction terms with the different impurity/Bose 
interaction strengths. The motion appears to be nearly that of a free particle, except the 
velocity of the particle is different from the free particle value po/M. The plot shows two 
curves with different particle/BF interaction strength, (3 = 0.3, and f3 = 1.0. The stronger 
interaction leads to a wavy motion of the particle. The oscillation of energy between the 
particle and BF can be attributed to the choice of initial conditions for BF, no = N, n\. = 
for k ^ 0, which are not in the equilibrium with the moving particle. This result clearly 
illustrates the microscopic superfluidity of such system. The particle cannot transfer energy 
to the Bose system, yet the interaction with Bose liquid causes the particle to have an 
effective mass M e ff = p /(dx(t)/dt). 

To better understand the dynamics observed in these calculations let us first consider the 
dynamics of a pure Bose system (Eqs. J9|10D with Vj = 0). Of course, the minimum of the 
Hamiltonian function is the stationary solution of equations of motion. Let us now consider 
a small deviation of one of the momentum coordinates from zero ^ and a* k ^ 0. Doing 
so one obtains coupled equations of motion for and a* k 

( k 2 n U k \ n Q , 
- ™ k = + —J a k + -U k a__ k , (12) 

/ k 2 n U k \ , n Q 

«o* = + ~ir) ak + 7T • (13) 

The simple diagonalization of these equations leads to the two normal modes that have the 
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frequency given by the Bogoliubov spectrum 

w-ft&W-Wf- <i4) 

Thus the "classical" frequency spectrum of the BF normal modes is the Bogoliubov's spec- 
trum. Note, that although we used the condensate number of particles, n = aga , in the 
above equations it is essentially the same as the total number of particles, n ~ N, due to 
the chosen initial conditions. 

A particle moving with a velocity v exerts the time dependent driving force on every of 
the normal modes of the BF given by 

L L 

This force is off-resonance with every normal mode of the BF if kv < e(k), i.e. the Landau 
criterion is satisfied. Nevertheless, this force results in the forced oscillations of the normal 
modes giving rise to rescaling of the particle velocity. The energy at given velocity is the 
energy of the particle and the total energy of all forced oscillations of the Bogoliubov's 
excitations. 

This classical picture is rather simple and should translate directly into the language 
of perturbation theory of quantum mechanics. However, the main difference between the 
Gross-Pitaevskii equations of motion and quantum mechanical solution turns out to be the 
ground state of the BF. The minimum of the classical Hamiltonian is characterized by all 
particles being in the k = mode. In the vicinity of this point the spectrum can only depend 
on the interaction with k = mode. In quantum mechanics this cannot be the case. The 
discrepancy can be thought to arise from non-commutation of the and operators or 
the zero point vibrations of Bogoliubov's excitations. This makes an anharmonic interaction 
term rather large even for a single quantum excitation of the BF. 

IV. EXCITATION SPECTRUM OF PURE BOSE FLUID 

In this section the full quantum mechanical problem of pure Bose fluid is solved by direct 
diagonalization of Hamiltonian Eq. 0. Since this Hamiltonian commutes with the operator 
of total momentum, it is sufficient to consider the blocks of the full Hamiltonian at every 
total momentum k. The Hamiltonian matrix is constructed in a given basis set and its blocks 
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at each k are diagonalized. The analysis reveals which one of the eigenstates (often, but 
not always, the lowest energy state) corresponds to the elementary excitation. Its energy is 
chosen to be the excitation energy E(k). The calculations are then repeated for all relevant 
values of k. 

Two schemes were implemented to perform such calculations. The first approach directly 
employs the space of occupation numbers as a basis set. Basis states are given as |{nfc}) We 
then truncate the highly excited states and retain only the states that have n > n min , where 
n min is a chosen parameter. The convergence of the results with n min is verified by changing 
this number. The implementation of this scheme leads to rather expensive calculations and 
the convergence was achieved only for weakly interacting systems (a ~ 0.1) or very small 
systems (on the order of 10 momentum modes total). 

The second approach involves a use of a significantly more efficient basis. Knowing that 
the normal modes of the classical dynamics are the Bogoliubov's excitations we develop 
a basis set based on the Bogoliubov's excitations rather than the occupation numbers of 
momentum wavefunctions. Let us use the main property of the Bose condensate, the fact 
that (no) is the macroscopic number (or simply large in the case of a finite system). 

The first step is to introduce new set of creation and annihilation operators 

a^a k , (16) 



m 

At = -±=a ai . (17) 
v n o 



These operators do not actually change the number of particles but rather promote particles 
from to A; momentum states. The commutation rules for these operators can be easily 
obtained to be 

[A k , At] = (1 + n k /n ) 5 k , k > ~ 5 k;k , . (18) 

which reduces to the regular Bosonic commutation in the limit of large no- We can then 
insert the identity operator / = aJa /no to rewrite the Hamiltonian via the new operators 
as 

* = ¥ + £ [(£ + ^ ) *t*>+ n 4r + ^4 <-> 

+ ^ J2 Ul ( A t A t-i A k + A t A k-iA^ + J2 U kl _ k3 AtAt 2 A k3 A k4 S kl+k2 _ k3 _ k4 . 
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Note that we effectively excluded the creation/annihilation of the particles with k = 0. The 
sums in the Hamiltonian go over all momenta k ^ and the Hamiltonian breaks into terms 
that are quadratic, cubic, and fourth-oder in these operators. The quadratic terms, that 
constitute the Bogoliubov's Hamiltonian, consist of the kinetic energy and the part of the 
interaction that is proportional to uq. The cubic terms are proportional to ^/uq and the 
quartic terms are of the order of unity 

The next step is to diagonalize the quadratic part of the Hamiltonian. Following Bogoli- 
ubov we introduce the new operators using relationships, 

A k = f l— J {B k + L k B + _ k ) , (20) 

V 1 ~ L k 

A + k = S^ {B+ + L k B_ k ) . (21) 

V 1 ~ L k 

The coefficient L k is chosen such that quadratic part becomes diagonal. As it is well known 
this procedure leads to the Hamiltonian, 

H = Y1 <k)B+B k + ({J3+ B k }) + ffW ({J3+ B k }) , (22) 



with the spectrum e(k) described by Eq. (14) and the coefficients 



1 ■ (23) 



n U k V 2M L 

We do not neglect any of the higher order terms which, as calculations show, have significant 
effect on the spectrum. The terms that are cubic and 4-th order in operators B k , B£ are rep- 
resented by and respectively. They can be obtained by substituting Eqs. ( 20|21 ) 



into Eq. (19). For the sake of space we do not explicitly write down those obvious but 
lengthly expressions. 



It readily verified that the ground state of the quadratic part of Hamiltonian Eq. (22) is 
given by 



\0) B = exp 



. k 



|0> , (24) 



where |0) is the fully condensed ground state of ideal Bose gas (n = N). This is the 
state with no Bogoliubov's excitations. The basis set is then generated by the action of 
Bogoliubov's creation operators onto this state. 

\n kl ,...,n kl ) = ^=(B+) |0) b . (25) 
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The numerical code generates a list of such basis states and sorts it according to their total 
momentum. A truncation of this basis at the maximum of 5 excitations is shown to produce 
converged results for a wide range of system parameters. A code then uses Hamiltonian 



Eq. (22) to build the Hamiltonian matrices for each total momentum k of the system. 



These matrices are then diagonalized resulting in a set of eigenenergies, Eq,E^" , ... , and 
corresponding eigenstates, l^o^)) )> ••• • 

At low total momentum k the lowest energy, Eq\ is the energy of the elementary ex- 
citation; however, at higher k this is not true, as the elementary excitation energy should 
approach k 2 /2m and become higher in energy than multiple excitations that combine to the 
total momentum k. Thus one needs a criterion to select an elementary excitation energy out 
of all eigenvalues. A simple strategy is to use the value of the overlap |(^ |S^"|0)b| 2 ; how- 
ever, for the strongly interacting BF all of such overlaps become very small as both ground 
and excited states deviate strongly from the Bogoliubov's states. A physics suggests another 
criterion that is used in this work. By examining the interaction terms in the Hamiltonians, 
Eq. {Eg, one realizes that excitations should be classified according to the overlap, 



I<#°|P*|4 0) >| 2 - (26) 

Indeed, this matrix element governs participation of an excited state in most physical pro- 
cesses. Note that the right hand side of this expression is precisely the state used by Bijl 
and Feynman to develop the theory of helium excitation spectrum [22]. 

The calculated spectrum is shown in Fig. 3. The parameters of the system are taken 



be: R = 3.0, N = 19, a = 0.5, and m — 1. The Bogoliubov's spectrum Eq. (14) as well 
as the k 2 /2m curve are shown for comparison. The deviation of the spectrum from the 
Bogoliubov's result is significant at low k. At small momenta (k < 5/3) the lowest energy 
state coincides with the elementary excitation. At intermediate momenta the single and 
double excitations become strongly mixed. This results in several states having significant 



overlap Eq. (26). Thus several points are shown in Fig. 3 all of which are the excited states 
of the system with momentum k but they cannot be classified as single or double excitations. 
Finally, at large momenta the elementary excitations approach the k 2 /2m limit of a single 
particle moving with the momentum k independently of BF. 

Despite the fact that the system is one-dimensional its excitation spectrum clearly shows 
linear rise at small momentum, thus this system should retain superfluid-like properties for 
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the purpose of impurity dynamics at zero temperature. The calculated slope of the spectrum 
at small k is significantly different than that of a Bogoliubov's spectrum, which indicates 
limited applicability of the classical mean-field treatment presented in previous Section. In 
the next Section we add an impurity to the BF and compute its energy spectrum using the 
results of the calculations presented above. 

V. ENERGY SPECTRUM OF IMPURITY 

It is natural to use the eigenstates of the BF obtained in the previous section in order 
to calculate the energy spectrum of an impurity. Using the momentum representation for 
an impurity particle we construct a basis set from the momentum states of an impurity and 
eigenstates of the BF found in the previous section. A basis function for a motion with total 
momentum p can be written as 

|$jg) = e %-*)*|^ fc )) . (27) 

Momentum, p, in this expression is a good quantum number; it is distributed between the 
impurity moving with momentum p — k and the excited state of BF with momentum k. We 
truncate the number of eigenstates of the BF that are used in the actual calculations. It 
is verified that the use of the 20 lowest states for each momentum k, (j = 1,2,... ,20) is 
sufficient; no improvement has been found for larger basis sets. 

In order to calculate the matrix of the full Hamiltonian the following matrices of the 
density operator are pre-computed and saved as the result of the pure BF calculations 
described in the last section, 

= (4 k) \P k -i\^) ■ (28) 

Using the Hamiltonian of the impurity in the form of Eq. (|6| the Hamiltonian matrices are 
computed as 

<!, = <*£ 4 I#6/ + H imp \^ d ) = + 4 kl) ) <W 2 <% + \v kl „ hlP ^ . (29) 

The energy spectrum of impurity E(p) is found by diagonalization of these matrices. The 
results are shown in Fig. 4. The parameters of the system are: R = 3.0, N = 19, a = (3 = 0.5, 
and M = m = 1. The k 2 /2M spectrum of the free particle is shown for comparison. The 
dotted line is the k 2 /2M e ff curve that is fitted to the first calculated point, p — 1/R. This 
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leads to M e // = 1.29, compared to M = 1. This line passes through the second calculated 
point with the 0.2% accuracy, after which the deviation of the calculated spectrum from the 
effective mass approximation becomes visible. 



VI. RELATION OF EXCITATION SPECTRA TO THE STRUCTURE OF THE 
GROUND STATE 

Computational results presented in previous sections clearly demonstrate the microscopic 
superfluidity in the finite model system. The excitation spectrum of an impurity consists of 
true eigenstates that are well separated from the rest of the excited states of the system at 
small k. Interestingly, the finite nature of the model does not have any effect on the physics 
except for the quantization of momentum. The initial several points of the impurity excita- 
tion spectrum are in nearly perfect agreement with the quadratic fit, i.e. the effective mass 
approximation. Our calculations show that this behavior is true regardless of the parameters 
of the system such as density of Bose particles and the strength of their interaction. These 
results appear to be a successful numerical experiment. While this experiment clearly shows 
microscopic superfluidity in a finite system, the results do not directly give insight into the 
underlying physics beyond the Bogoliubov's approximation. It is clear, however, that the 
latter gives the poor estimate of the effective mass as well as poor physical description of 
the BF. In this section we present a microscopic theory that explains the nature of the BF 
spectrum as well as an impurity energy based on the structural properties of the ground 
state. 

A successful theory of the superfluid helium excitation spectrum has been build on the 
assumption that the properties of the ground state govern the energy of the excitations. In 
particular, Feynman argued [23] that given the unknown ground state wavefunction \ipo), 
the most natural wavefunction of the excited state with momentum k is given by 



He then used the properties of the ground state to evaluate the average energy of the BF in 
such state. The first order approximation to the energy spectrum becomes 



(30) 



E(k) 



k 2 1 
2^S(fc) ' 



(31) 
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where S(k) is the ground state structure factor defined by 

S(k) = -^(Mp-kpk\%) ■ (32) 

The structure factor can be related directly to the radial distribution function (RDF), g(r), 
of the BF as 

g(r) = l-J2e lkr (S(k)-l). (33) 

k 

This theory can easily be verified using the results of the calculation of section IV. The 



RDF, Eq. (33), calculated using the computed ground state of BF is shown in Fig. 5. It 
is interesting to note that the for the particular choice of system parameters the repul- 
sion of the Bose particles is clearly observed in the RDF as a well at r = 0; however, the 
probability of two particles occupying the same volume is significant, in fact the density 
of BF at zero distance from a particle is (p(0)/(p) = 0.37). As we found out, the calcu- 
lations become prohibitively expensive when the displacement of BF by an atom becomes 
complete, i.e. p(0) — > 0, because the states with large number of Bogoliubov's excitations 
contribute significantly to the ground state. The Bijl-Feynman spectrum is shown in Fig. 6 
and compared to the calculated one. As in the case of the superfluid helium, the spectrum 
is in excellent agreement at low values of k. At intermediate momenta the Bijl-Feynman 
expression overestimates the excitation energy. The reason of this effect is clear and can be 



understood from numerical results: the wavefunction Eq. (30) becomes strongly mixed with 
double excitations thus lowering the real excitation energy. 

Similar theory has been developed for the excitation spectrum of impurity [25J. Let us 
assume that |$o) is the ground state of the combined impurity-BF system. Wavefunction 
that represents motion of impurity independent of BF is given by 

|$<W) = e ^|$ ) . (34) 

As expected, the variational energy of this wavefunction gives no correction to the mass of 
the particle resulting in the energy, 

The next step is to incorporate the possible excitations of the BF into the motion of impurity. 
While the exact wavefunctions of BF excited states are not known it is easy to incorporate 
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the Feynman-like excitations by using the states 



1*1") 



1 



(36) 



^/NS{l) 

The pre-factor (iVS^Z)) -1 / 2 can be shown to normalize this wavefunction. The state of 



Eq. (36) has the total momentum k, which results from a particle moving with momentum 
k — I and BF state of momentum I. These states provide a very natural basis set for the 
moving particle giving the expression for the moving particle as 



(37) 



Such states are not orthogonal and their overlap can be evaluated; however, it is does not 
enter the perturbation expression derived below. An important property of these states is 
that the relevant Hamiltonian matrix elements can be easily evaluated and expressed through 
the structural properties of the ground state. Using the method developed by Feynman one 
obtains 



(0\ 



and 



(0)v 



(fc-Z) 2 I 2 
2M + 2mS{l) 

kl G{1) 



Here we introduced the impurity/BF structure factor 



G(k) = ($ |e- fe p fc |$o) 



(38) 
(39) 

(40) 



Using Eqs. ( 38|39 ) the second order perturbation correction to the energy Eq. (35 ) is obtained 



E 



(2) 



(Jd_\ 2 G(Q 2 
\2M) NS(l) 



k 2 
I 2M 



(k-iy 

2M 



V 2 



(41) 



2mS(l) 



At small k, the main term of the denominator Eq. (41) is the energy of the BF excitation 
l 2 /[2mS(l)], thus this equation can be reduced to 



k 2 m 1 



E 



M 2 N \S(l) J 
which leads to the effective mass of impurity given by 



(42) 



M, 



eff 



M 



m 1 . ^ (G{1)\ 



(43) 
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Once again, this theory can be easily tested using our numerical results. The impurity/BF 
structure factor G(k) is essentially a Fourier transform of the BF density distribution around 
the particle. Both functions are obtained by the analysis of the ground state obtained in 
our calculations. The density distribution of Bose particles around an impurity, given by 



p(r) 



\\N-^G{ky k ^j , (44) 



is shown in Fig. 7. This function is similar to the RDF of BF, except is not normalized to 
unity at large distances but instead represents the actual density of the BF. In our particular 
case the density of Bose particles is increased in the region away from impurity compared 
to the unit value because of the actual displacement of particles in the small system. The 



energy spectrum predicted by Eqs. (41) is compared to the calculated spectrum in Fig. 8. 
The result shows nearly perfect numerical agreement with the first two points of the spectrum 
and starts to deviate at larger k. The main source of the deviation can be identified as the 



interaction with double excitations not considered by Eq. (37). The value of effective mass 



predicted by Eq. (43), M e ff/M = 1.27, is very similar to the value 1.29 computed directly. 
It is interesting to note, that the effective mass cannot be understood simplistically as the 
mass of BF displaced by the interaction with impurity particle. As our calculations show, 
M displaced ~ 2, while M eff -M = 0.29. 



The perturbation theory expression, Eq. (41 ) is equivalent to the one-phonon intermediate 



perturbation (OIP) result obtained in ref. [25]. Although this equation is not as accurate 
for a more realistic system of 3 He in 4 He, our calculations show that it captures the essential 
physics of impurity motion. For our system this result gives an accurate description of 
impurity excitation spectrum at small k. The most important feature of this theory is 
that excitation energies are expressed entirely through the structure of the ground state. 
The latter can be computed for a realistic system such as liquid He using the wealth of 
the imaginary time path integral techniques. In other words, the dynamical problem is 
reduced to a statistical one. A rather complete effort in this direction has been recently 
undertaken by Zillich and Whaley [T5]. Their work combines the variational CBF theory 
with the diffusion Monte Carlo simulations to compute moments of inertia of molecules in 
helium droplets. Our work shows quantitative success of much simpler approximations based 
entirely on the density distributions. This success raises hope that simpler methods based 
on the perturbation expansions could be used to solve similar problems. 
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VII. CONCLUSIONS 



In this work we present a computational study of a finite one-dimensional system that 
consists of a BF and an additional impurity particle confined to a ring. None of the inter- 
actions are taken to to be small yet the full quantum calculations are converged for this 
system. The effects of the strong particle repulsion is clearly observed in the structure of the 
ground state. It is shown that this system exhibits a property of microscopic superfluidity, 
i.e. there is a brunch of energy states that corresponds to the motion of an impurity particle. 
These are true energy eigenstates of the system that are well separated from the continuum 
of quantum states of the BF. The calculations suggest a most natural analytical description 
of impurity motion. In particular, the excited states of the impurity can be obtained from 
the ground state of the impurity-BF system. It is shown that the excitation spectrum can 
be accurately predicted from the structure of the ground state. 

We believe that the calculations presented in this work provide a valuable illustration 
to the analytical theory of superfluidity developed over the last sixty years. The classical 
dynamics nearly ideally illustrates the work of Bogoliubov exemplifying the main approx- 
imation inbuilt in his treatment. The dynamics of this system clearly shows the validity 
of the Landau criterion and shows that it is applicable to finite systems. Finally, the full 
quantum calculations illustrate in detail the microscopic theory of superfluidity developed 
by Feynman and his successors. It is our hope that a complete understanding of this model 
system will help us to combine the modern state of the art computational methods with the 
wealth of analytical theory developed decades ago. 
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Figures 




FIG. 1: A finite one-dimensional model of Bose fluid (BF). Bose particles are represented by an 
empty circles. An additional impurity particle is represented by the filled circle. The calculation 
of the BF and the impurity excitation spectra is the subject of this paper. 
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FIG. 2: . Classical trajectories of impurity interacting with BF obtained by solving Eqs. ([7 10). 
The dashed line represents the free-particle dynamics. Top plot shows the motion of impurity 
interacting with an ideal Bose Gas (a = 0). The bottom plot shows the motion impurity through 
the interacting BF (a = 0.5); two curves are the trajectories that correspond to two different 
impurity/BF interaction parameters ((3 = 0.3 and (3 = 1.0). 
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FIG. 3: Crosses show the calculated excitation spectrum, E(k), of pure BF obtained by full 
quantum mechanical solution of Hamiltonian Eq. ([T|. The Bogoliubov's excitation spectrum, 



Eq. (14), (dashed line) and k 2 /2m curve (dotted line) are shown for comparison. 
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FIG. 4: Crosses show the calculated excitation spectrum of an impurity particle, E(k). The 
solid line is the k 2 /2M curve, with M = 1. The dotted line is the k 2 /2M e ff fit to the spectrum. 
M e ff = 1.29 is obtained using the lowest point on the spectrum. 
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FIG. 5: Radial distribution function, g(r) of the BF computed using the ground state calculated 
in section IV. 
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FIG. 6: The Bijl-Feynman spectrum, E(k) = k 2 /[2mS(k)], (solid line) is compared to the spectrum 
calculated in section IV. S(k) is the structure factor of the ground state computed in section IV. 
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FIG. 7: Average density of BF as a function of distance from impurity p{r) resulting from the 
ground state of impurity/BF system computed in section V. 
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FIG. 8: Impurity excitation spectrum calculated using Eq. (41 ) (circles) compared to the spectrum 



calculated in section V (crosses). Dotted line is the effective mass approximation. 



26 



